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Using renormalization group methods we evaluate the thermal strangeness and charm chemical equilibrium 
relaxation times in the deconfined quark-gluon plasma. We present a reaction model and evaluate the total 
production rate of strangeness in fixed target Pb— Pb collisions at 10—300 A GeV. We discuss the relevance of our 
results to the diagnosis and understanding of the properties of the deconfined state. 



1 Introduction 

In high energy nuclear collisions^! we are interested 
to create and study the deconfined quark-gluon 
phase, believed to consist in the limit of very high 
energy densities of a liquid of quarks and gluons 
interacting perturbatively. We are motivated by 
the desire to recreate in the laboratory conditions 
akin to those prevailing in the first moments of 
the Early Universe, and in the recognition that we 
can study the properties of the excited, 'melted', 
vacuum state of strong interactions. 

Since in the collision of large nuclei the highhj 
dense state is formed for a rather short timcq 
one of the major challenges has been to identify 
suitable physical observables of deconfinement. A 
number of possible experimental signatures of the 
formation and properties of quark-gluon plasma 
(QGP) have been studied. This report addresses 
our recent advances in evaluating the relevance pi 
strangefju and charm quark flavor in this questo. 
Other major probes include the phenomenon of 
Jfl! suppression B, photons and dileptonstffl. 

Strangeness is a very interesting diagnostic 
tool of dense hadronic matterB: 

1) particles containing strangeness are found 
more abundantly in relativistic nuclear collisions 
than it could be expected based on simple scaling 
of p—p reactions; 

2) all strange hadrons have to be made in in- 
elastic reactions, while light u, d quarks are also 
brought into the reaction by the colliding nuclei; 

3) because there are many different strange 



particles, we have a very rich field of observables 
with which it is possible to explore diverse prop- 
erties of the source; 

4) theoretical calculations suggest that glue- 
glue collisions in the QGP provide a sufficiently 
fast and thus by far a unique mechanism leading 
to an explanation of strangeness enhancement. 

We begin by recalling these mechanisms of 
strangeness production. We then obtain in section 
|^ the magnitude of running coupling constant and 
quark masses. Using these results we evaluate the 
thermal relaxation times of strangeness and charm 
as function of temperature in section |[ We dis- 
cuss the relevance of the flavor observable of QGP 
and describe how these relaxation times allow to 
compute the hadronic particle yields in section [|. 
We close with a few general remarks about the 
relevance of our results. 



2 Strangeness Production in QGP 

We will now show how to evaluate using twn_par- 
ticle collision processes flavor productionErtlj in 
thermal QGP. Ultimately, we will employ run- 
ning QCD parameters a s (/i) and rrii(fj,), i =s,c. 
While theory and experiment constrain now suffi- 
ciently the coupling strength a s , considerable un- 
certainty still remains in particular in regard of 
strange quark mass scale, as well as systematic 
uncertainty related to applications of QCD to soft 
(less than 1 GeV) processes. 

The generic angle averaged two particle cross 
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section for (heavy) flavor production processes 
9 + 9 — ► / + / and q + q -> / + / ' , are: 



'gg->ff 



f(s) 



27ra 
~37 



2 r 



f 



4mf 



(1) 



tanh _1 W(s) - 



31mf 



W(s) 



(s) 



87T0| 

27s 



2m\ 



W(s) 



(2) 



where W(s) = \/l — Arrif/s , and both the QCD 
coupling constant a s and flavor quark mass mi 
will be in this work the running QCD parameters. 
In this way a large number of even-ce s diagrams 
contributing to flavor production is accounted for. 

What remains unaccounted for is another 
class of processes in which at least one additional 
gluon is present. In particular processes allowing 
the production of an additional soft gluon in the 
final state remains unaccounted for today. Lead- 
ing diagrams contain odd powers of a s and their 
generic cross section is in general infrared diver- 
gent, requiring a cut-off which for processes oc- 
curring in matter is provided by the interactions 
(dressing) with other particles present. The pro- 
cess in which a massive 'gluon', that is a quasi- 
particle with quantum numbers of a gluon, decays 
into a strange quark pair, is partially included in 
the resummation that we accomplish in the pres- 
ent work. At the present time we do not see a 
systematic way to incorporate any residue of this 
and other effects, originating in matter surround- 
ing the microscopic processes, as work leading to 
understanding of renormalization group equations 
in matter (that is at finite temperature and/or 
chemical potential) is still in progresscil 

3 Running QCD Parameters 

To determine the two QCD parameters required, 
we will use the renormalization group functions /3 
and 7, n : 

doi dfrt 
^~q~ = 0( a s(v)) , = -m7 m (a s (^)) . (3) 

For our present study we will use the perturbative 
power expansion in a s : 



/?P crt = a 2 s [ b + ha s + b 2 a\ 
7 m ert = a s [ c + c\a s + C2a^ 
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For the SU(3)-gauge theory with m fermions the 
first two terms (two 'loop' order) are renormaliza- 
tion scheme independent, and we include in our 
calculations the three 'loop' term as well, which is 
renormalization_scheme dependent, evaluated in 
the MS-schemeEl We have: 
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The number rif of fermions that can be excited, 
depends on the energy scale \i. We have imple- 
mented this using the exact phase space form ap- 
propriate for the terms linear in nf. 
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with m s = 0.16 GeV, m c = 1.5 GeV, m b = 4.8 GeV. 
We checked that there is very minimal impact of 
the running of the masses in Eq. (7) on the fi- 
nal result, and will therefore not introduce that 
'feed-back' effect into our current discussion. The 
largest effect on our solutions comes from the bot- 
tom mass, since any error made at about 5 GeV is 
amplified most. However, we find that this results 
in a scarcely visible change even when the mass is 
changed by 10% and thus one can conclude that 
the exact values of the masses and the nature of 
flavor threshold is at present of minor importance 
in our study. 

We show the result of numerical integration 
for a s in the top portion of Fig. 0. First equation 
in (||) is numerically integrated beginning with an 
initial value of a s (Mz). We use in this report the 
August 1996 World averageB a s {M z ) = 0.118 
for which the estimated error is ± 0.003 . This 
value is sufficiently precise to eliminate most of the 
uncertainty that has befallen much of our earlier 
studicsO'Ej. In addition, the thin solid lines present 
results for a s {Mz) = 0.115 till recently the pre- 
ferred result in some analysis, especially those at 
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running mass equation is linear in m, it is possible 
to determine the universal quark mass scale factor 
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Figure 1: a B (fi) (top section); the equivalent parameter Ao 
(middle section) and m T (fj.) = m(fM)/m(Mz) (bottom sec- 
tion) as function of energy scale fi. Initial value a B (Mz) = 
0.118 (thick solid lines) and a s (M z ) = 0.115 (thin solid 
lines). In lower section the dots indicate the strangen- 
ess pair production thresholds for m B (Mz) = 90 MeV, 
while crosses indicate charm pair production thresholds for 
m c (M z ) = 700 MeV. 



lower energy scale. As seen in Fig. |T|, the variation 
of a s with the energy scale is substantial, and in 
particular we note the rapid change at and below 
fj, = 1 GeV, where the strange quark flavor forma- 
tion occurs in hot QGP phase formed in present 
day experiments at 160-200 A GeV (SPS-CERN). 
Clearly, use of constant value of a s is hardly justi- 
fied, and the first order approximation often used: 
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leads to a strongly scale dependent Aq(^) shown 
in the middle section of Fig.[l]. Thus it also cannot 
be used in the evaluation of thermal strangeness 
and charm production. 

With a s (/i) from the solutions described 
above, we integrate the running of the quark 
masses, the second equation in p). Because the 



m r = m(ji)/m{fio) . 
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Since a s refers to the scale of fio = Mz, it is a 
convenient reference point also for quark masses. 
As seen in the bottom portion of Fig. |l|, the change 
in the quark mass factor is highly relevant, since 
it is driven by the rapidly changing a s near to 
fi ~ 1 GeV. For each of the different functional 
dependences a s ([i) we obtain a different function 
m I . The significance of the running of the charmed 
quark mass cannot be stressed enough, especially 
for thermal charm production occurring in fore- 
seeable future experiments well below threshold, 
which amplifies the importance of exact value of 
m c . 

Given these results, we find that for a s = 
0.118 and m s (M z ) = 90 ± 18 MeV a low energy 
strange quark mass m s (l GeV) ~ 200 ±40 MeV, in 
the middle of the standard range 100 < m s (l GeV) 
< 300 MeV. Similarly we consider m c (Mz) 
7i 10 ± 50 MeV, for which value we find the low 
energy mass m c (l GeV) ~ 1550 ± 110 MeV, at the 
upper (conservative for particle production yield) 
end of the standard range 1 < m c (lGeV) < 1.6 
MeV. There is another (nonperturbative) effect of 
mass running, related to the mass at threshold for 
pair production m* h , i = s, c, arising from the so- 
lution of: 
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This effect stabilizes strangeness production cross 
section in the infrared: below y/s = 1 GeV 
the strange quark mass increases rapidly and the 
threshold mass is considerably greater than m s (l 
GeV). We obtain the threshold values 2m* h = 611 
MeV for a s (M z ) = 0.118 and 2m* h = 566 MeV for 
a s (Mz) = 0.115. Both values are indicated by the 
black dots in Fig.[l]. For charm, the running mass 
effect plays differently: since the mass of charmed 
quarks is listed in tables for fi = 1 GeV, but the 
value of the mass is above 1 GeV, the production 
threshold mass is smaller than expected (i.e., listed 
value). For m c (Mz) = 700 MeV the production 
threshold is found at ~ 2m* h ~ 2.3 GeV rather 
than 3.1 GeV that would have been expected for 
the m c (l GeV). This reduction in threshold en- 
hances thermal production of charm, especially so 
at low temperatures. 
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Figure 2: QGP strangeness relaxation time, for a B (Mz) = 
0.118, (thick line) and = 0.115 (thin line); m s (M z ) = 
90 MeV. Hatched areas: effect of variation of strange quark 
mass by 20%. Dotted: comparison results for fixed a s = 0.5 
and m s = 200 MeV. 



4 Strangeness and Charm Thermal Relax- 
ation Times 

The thermal average of the cross section is the in- 
variant production rate per unit time and volume: 
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The dots indicate that other mechanisms may 
contribute to strangeness production. The par- 
ticle distributions are in our case thermal 
Bose/Fermi functions (for fermions with A q = 
1.5), and g q = 6, g & = 16. For strangeness pro- 
duction rif = 2, and for charm production rif = 3 . 
^From the invariant rate we obtain the strangen- 
ess relaxation time r s shown in Fig.|^, as function 
of temperature: 
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Note that here unaccounted for processes, such as 
the above mentioned odd-order in a s would add 
to the production rate incoherently, since they 
can be distinguished by the presence of incom- 
ing/outgoing gluons. Thus the current calculation 
offers an upper limit on the actual relaxation time, 




t- 10 2 - 



_'()() 



ooo 



100 



000 



600 



T [MeV] 



Figure 3: Solid lines: thermal charm relaxation constant 
in QGP, calculated for running a s {M z ) = 0.115; 0.118, (in- 
distinguishable), m c (Mz) = 700 MeV. Lower dotted line: 
for fixed m c = 1.1 GeV, a a = 0.35; upper doted line: for 
fixed m c = 1.5 GeV, ct B = 0.4. Hatched area: effect of 
variation m c (M z ) = 700 ± 50 MeV 



which may still be smaller. In any case, the present 
result suffices to confirm that strangeness will be 
very near to chemical equilibrium in QGP formed 
in collisions of large nuclei. 

We show in Fig.|| also the impact of a 20% 
uncertainty in m s {Mz), indicated by the hatched 
areas. This uncertainty is today much larger com- 
pared to the uncertainty that arises from the re- 
cently improved precision of the strong coupling 
constant deterrcunationlla. We note that the cal- 
culations madea at fixed values a s = 0.5 and 
m s = 200 MeV (dotted line in Fig.|) are well 
within the band of values related to the uncer- 
tainty in the strange quark mass. 

Since charm is somewhat more massive com- 
pared to strangeness, there is still less uncertainty 
arising in the extrapolation of the coupling con- 
stant. Also the systematic uncertainty related to 
the soft gluons (odd-a s ) terms are smaller, and 
thus the relaxation times r c we show in Fig.|^ 
are considerably better defined compared to r s . 
There is also less relative uncertainty in the value 
of charm mass. We also show in Fig. | (dotted 
lines) the fixed m c , a s results with parameters se- 
lected to border high and low T limits of the re- 
sults presented. It is difficult to find a good com- 
parative behavior of r c using just one set of m c and 
a s . This may be attributed to the importance of 
the mass of the charmed quarks, considering that 
the threshold for charm production is well above 
the average thermal collision energy, which results 
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in emphasis of the effect of running charm mass. 
In the high T-limit the choice (upper doted line 
in Fig.||) m c = 1.5 GeV, a s = 0.4 is appropri- 
ate, while to follow the result at small T (lower 
doted line in Fig.|^) we take a much smaller mass 
m c = 1.1 GeV, a s = 0.35. 

We recall that the equilibrium distribution is 
result of Boltzmann equation description of two 
body collisions. Thus the mass arising in the equi- 
librium density in Eq. (|l^) is to be taken at 
the energy scale of the average two parton col- 
lision. We adopt for this purpose a fixed value 
fh s — 200 MeV, and observe that in the range of 
temperatures here considered the precise value of 
the mass is insignificant, since the quark density 
is primarily governed by the T 3 term in this limit, 
with finite mass correction being 0(10%). The 
situation is less clear for charm relaxation, since 
the running of the mass should have a significant 
impact. Short of more complete kinetic treatment, 
we used m c ~ 1.5 GeV in order to establish the 
reference density in Eq. (|l^) . 

5 QGP Flavor Observable 

We will indicate in this section how the study of 
flavor production impacts our understanding and 
diagnosis of the deconfined QGP phase. We recall 
first that there are two generic flavor observable 
which we can study analyzing experimental data: 

• yield of strangeness/charm: 

once produced in hot early QGP phase, strangen- 
ess/charm is not reannihilated in the evolution of 
the deconfined state towards freeze-out, and thus 
the flavor yield is characteristic of the initial, most 
extreme conditions; 

• phase space occupancy 7 S , C : 

impacts distribution of flavor among final state 
particle abundances. 

Given that the thermal equilibrium is estab- 
lished within a considerably shorter time scale 
than the (absolute) heavy flavor chemical equili- 
bration, we can characterize the equilibration of 
the phase space occupancy by an average over the 
momentum distribution: 

J d 3 pd 3 xni(p,x;t) . 
J d?p cfo nf° (p, x) ' 

The chemical equilibrium density is indicated by 
upper-script 'oo'. When several carriers of the 



flavor are present, as is the case in the confined 
phase, Hi is understood to comprise a weighted 
sum. 

In order to be able to compute the produc- 
tion and evolution of strangeness and charm fla- 
vor a more specific picture of the temporal evo- 
lution of dense matter is needed. Here, we will 
address specifically strangeness production in col- 
lisions at CERN-SPS, up to 200 A GeV per nu- 
cleoli. We use a simple, qualitative description, 
simplified by the assumption that the properties 
of the hot, dense matter are constant across the 
entire volume (fireball model). We consider ra- 
dial expansion to be the dominant factor for the 
evolution of the fireball properties such as tem- 
perature/energy density and lifetime of the QGP 
phase. Within the fireball model the expansion 
dynamics follows from two assumptions: 

• the (radial) expansion is entropy conserving, 
thus the volume and temperature satisfy: 

V ■ T 3 = Const. (14) 

• the surface flow velocity is given by the sound 
velocity in a relativistic gas 

v f = l/V3. (15) 

This leads to the explicit forms for the radius 
of the fireball and its average temperature: 

1 7] n 

R = R m H =(i — tin), T = r= . 

V3 K h l + (t-t in )/V3R in 

The initial conditions for Pb-Pb: 

A q = 1.6, t ln = lfm/c, = 320 MeV; with 

R in = 4.5 fm for n = 0.5; R m = 5.2 fm for n = 0.75, 

and for S-Pb/W: 

A q = 1.5, i in = lfm/c, ^ = 280 MeV; with 
Rin = 3.3 fm for 77 = 0.35; R ni = 3.7 fm for w = 0.5, 
have been determined such that the energy per 
baryon is given by energy and baryon flow, and 
the total baryon number is r)(Ai +A%), as stopped 
in the interaction region. The radius shown above 
are for zero impact parameter. For this, equations 
of state of the QGP are needed, and we have em- 
ployed our modeled in which the perturbative cor- 
rection to the number of degrees of freedom were 
incorporated along with thermal particle masses. 

In the fireball in every volume element we 
have: 

n s (p;t)=%n^(p;T^ s ). (16) 
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Table 1: 7 S and N B /B in S— W at 200 A GeV and Pb— Pb at 158 A GeV for different stopping values of baryonic number 
and energy rjg = rjg ; computed for strange quark mass m s (lGeV) = 200 ± 40 MeV, a B (Mz) = 0.118 . 





S-W at 200 A GeV 


Pb-Pb at 158 A GeV 


m = ve 


0.35 


0.5 


0.5 


0.75 


Is 


0.53 ±0.14 


0.65 ±0.15 


0.69 ±0.15 


0.76 ±0.16 


N s /B 


0.67 ±0.16 


0.70 ±0.16 


0.954 ±0.20 


0.950 ±0.20 



In this limit and allowing for the detailed balance 
reactions, thus re-annihilation of flavor, the yield 
is obtained from the equation: 



dN s (t) 
dt 



= V(t)A s [l- 7s 2 (i)] 



(17) 



Allowing for dilution-of the phase space density in 
expansion, we deriveu from Eq. fll7 ) an equations 
describing the change in %(t): 



dt 



Tm s 



rp 2 
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Here K2 is a Bessel function and x = m s /T. Note 
that even when 1 — 7^ < 1 we still can have a pos- 
itive derivative of 7 S , since the first term on the 
right hand side of Eq. ( |l8| ) is always positive, both 
T and d/ dx(x 2 k2) being always negative. This 
shows that dilution due to expansion effects in 
principle can make the value of 7 S rise above unity. 

Given the relaxation constant r s (T(t)), these 
equations can be integrated numerically, and we 
can obtain for the two currently explored experi- 
mental systems the values of the two observables, 
7 S and N s /B), which are given in table 1. 

For S-Ag collisions at 200 A GeV a recent 
evaluation of the specific strangeness yield leads 
to N s /B\ cxp = 0.86^ 0.14 (see table 4 of RefB). 
Our earlier analysiso of the WA85 data yields 7 S = 
0.75±0.1 for S-W interactions. Both these results 
are in good agreement with the theoretical result 
shown in the table, favoring the 50% stopping case 
for S-W. The analysis of the experimental Pb-Pb 
data is in progress. 

As we can see in table 1, there is a consid- 
erable uncertainty due to the unknown mass of 
strange quarks. On the other hand, inspecting 
the yield of strange quarks per baryon N s /B there 



seem to be very little dependence on the stopping 
fractions. This insensitivity to the reaction mech- 
anism coupled to a visible sensitivity to strange 
quark mass suggests that additional insight may 
be ultimately gained about the strange quark mass 
from the study of strangeness enhancement in rel- 
ativistic heavy ion collisions. Thus we show in 
Fig.|| as function of m s (Mz) for the two different 
stopping fractions the resulting strangeness yield 
per baryon, N$/B. The range of strange quark 
mass shown corresponds to the allowed range 45 < 
m s (Mz) < 135 MeV (100<m B (l GeV) < 300 MeV). 
The known systematic and statistical error is in- 
dicated by the divergence of the different curves, 
and is particularly small for N s /B. 

Before we can use these results to obtain a 
more reliable estimate of strange quark mass, we 
will need to understand other sources of system- 
atic uncertainty, e.g., those associated with un- 
known mechanisms of strangeness production and 
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Figure 4: Phase space occupancy 7 S and yield of strange 
quarks per baryon N B /B as function of strange quark mass 
m B (Mz) for Pb-Pb collision system at 158 A GeV, for the 
two stopping fractions and a B (Mz) = 0.118 ± 0.003 : thick 
solid lines, 77 = 75% , with long dashed lines indicating the 
a B (Mz) uncertainty, thick dashed lines for rj = 50% , with 
dotted lines indicating the a B (Mz) uncertainty. 
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E lob [A GeV] E, ob [A GeV] 

Figure 5: Phase space occupancy f B and yield of strange 
quarks per baryon N B /B as function of laboratory collision 
energy in Pb— Pb interactions for m B = 200 ± 40 MeV and 
two different stopping fractions: r\ = 50% (dashed and two 
doted lines) and 77 = 75% (solid and two long dashed lines). 



E ub [A GeV] E ub [A GeV] 

Figure 6: Relative yields of antibaryons (left) and bary- 
ons (right) as function of heavy ion collision energy i?Lab- 
Abundance set to a common value (unity) at i?Lab = 158 
A GeV. Solid lines nucleons, long dashed A + E°, short 
dashed S, dotted f2. 



to improve the model of QGP we are employing. 
However, these results are so strongly dependent 
on m s and so little on other quantities, that we 
can be optimistic that one day strange and charm 
flavor mass may be obtained from the flavor yields 
seen in nuclear collisions. 

Assuming that the model we proposed has 
been effectively tested at 200 A GeV S-W/Pb 
collisions, we compute the strangeness yield and 
phase space occupancy as function of energy. The 
results are shown for the case of Pb-Pb collisions. 
We show these results in Fig. || as function of lab- 
oratory collision energy, stressing that if at low 
energies the QGP phase is not encountered, we 
would expect to see a drop in strangeness yield 
beyond the expectations here presented. 

These results allow us to evaluate the strange 
(anti)baryon yields from QGP as function of col- 
lision energy. We note that at fixed m± the 
medium dependent factor controlling the abun- 
dance of hadrons emerging from the surface of the 
deconfined region is related to the chemical condi- 
tions in the source, and for strange quarks, there 
is also the occupancy factor 7 S to be considered: 

n h \ m± =e- m ^ /T Y[ lk X k . (19) 

The strange quark fugacity is in deconfined phase 
unity, while the light quark fugacity evolution 
with energy_pf colliding ions follows from our ear- 
lier studiesEj. In Fig.^. we have normalized all 
yields at i?Lab = 158 A GeV. Remarkably, all an- 
tibaryon yields (left hand side of Fig.||) cluster 



together (solid lines: (anti)nucleons, long dashed: 
(anti)hyperons, short-dashed: (anti)cascades, and 
dotted: (anti)omegas), thus as long as the QGP 
phase is formed, ratios of rare multistrange an- 
tibaryons should not change significantly while the 
collision energy is reduced, until the QGP forma- 
tion is disrupted. It should be noted that the yield 
of fl remains appreciable, all the way even at very 
small energies — this is the case as long as these 
particles are produced by the deconfined phase, 
rather than in individual hadronic interactions. 
For baryons (right hand side of Fig. ^) there is con- 
siderable differentiation of the yield behavior: the 
reference yield of nucleons rises slightly to compen- 
sate for the drop in the yield of other strange bary- 
ons, which decreases substantially with energy. 

6 Final Remarks 

In conclusion, we can once more remind our- 
selves about the two generic strangeness observ- 
ablcs. The relative total abundance of strangeness 
is most related to the initial condition, the 'hot- 
ter' the initial state is, the greater the produc- 
tion rate, and thus the final state relative yield, 
to be measured with respect to baryon number or 
global particle multiplicity (entropy). The phase 
space occupancy of strangeness 7 S depends aside 
of the initial production rate, on the final state 
dilution characterized by dynamics of the expan- 
sion and the freeze-out temperature. Our results 
show that we can use certain features of strange 
particle production to see the formation of the de- 
confined state and to study some QCD properties 
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and parameters. Our results suggest that already 
at present energies deconfinement is attained, and 
we have explored a number of features as function 
of collision energy in order to see if a more system- 
atic study is capable to confirm this conclusion. 

Using QCD renormalization group methods 
we have studied the s and c flavor chemical equi- 
librium relaxation times. We have shown that the 
newly measured QCD coupling constant comprises 
sufficiently small uncertainty to allow precise eval- 
uation of strangeness production at and below 1 
GeV energy scale. Our study has further proven 
that it is essential to incorporate in the evaluation 
of flavor production rates both running coupling 
constant and running mass. 

We find that running of the QCD parameters 
is of major significance, since, e.g., the effective 
charm production mass is considerably reduced, 
seen on the scale of available thermal energies. We 
found considerable enhancement of charm produc- 
tion for temperatures applicable at SPS collision 
energy, compared to fixed mass results. While 
charm experiences at low temperature T ~ 200 
McV a 100 times slower approach to chemical equi- 
librium compared to strangeness, for temperatures 
of about 500 MeV, as may apply to the conditions 
generated at LHC or perhaps even RHIC collider, 
t c — ► 30 fm, which is within factor two of the ex- 
pected maximum lifespan of the deconfined state. 
Thus our calculations suggest that there will be a 
significant abundance of thermal charm in nuclear 
interactions at RHIC/LHC. In consequence, open 
charm should play a similar role in the diagnosis of 
the 'hot' T ~ 500 MeV deconfined state as strang- 
eness is playing today for the 'cold' T ~ 250 McV 
case, and charm equilibrium appears within reach 
of the extreme conditions possibly arising at LHC. 

Our here presented results imply that in key 
features the strange particle production results ob- 
tained at 160-200 A GeV, are consistent with the 
QGP formation hypothesis. However, in order 
to ascertain the possibility that indeed the QGP 
phase is already formed today, a more system- 
atic experimental exploration as function of colli- 
sion energy of the different observable is required, 
for which purpose we also have explored the colli- 
sion energy dependence of the most characteristic 
strange particle features expected from the QGP 
phase. 
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